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Abstract. 

A simultaneous NLO QCD fit of the parton density functions (PDFs) of tfie proton and of the strong coupling as{Mz) is 
presented. The analysis is based on combined HI and ZEUS inclusive DIS cross section measurements together with inclusive 
jet cross section measurements provided by both HI and ZEUS collaborations. The inclusion of the jet data significantly 
reduces the correlation between the gluon PDF and the strong coupling, improving the precision of the gluon PDF and 
providing an accurate determination of as{Mz) using DIS data from HERA only. 
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INTRODUCTION 

The Rutherford centennial is celebrating Rutherford's discovery of the atomic nucleus in 1911 by scattering alpha- 
particles on a gold foil. Since then scattering experiments have come a long way, improving the resolution from about 
6 times the proton radius (Rp^i I fm) in 1911 to about 1/1 000 -Rpas achieved by the HERA experiments [ 1 ] at the so 
far only ep collider Here, I want to describe the cuiTent achievements in the extraction of the parton density functions 
(PDFs) of the proton and of the strong QCD coupling from inclusive and jet measurements in DIS at HERA. 

Inclusive neutral current (NC) and charged current (CC) measurements, in leading order (LO) represented by the 
quark-parton model, are sensitive to the valence and sea quarks of the proton. Only at next-to-leading order (NLO) and 
beyond do they become sensitive to the gluon density (scaling violations) and the strong coupling, which are however 
strongly correlated. Inclusive jets in DIS, when measured for example in the Breit frame, are sensitive to the gluon 
density and the strong coupling already in LO, via the boson-gluon fusion process. This process dominates for the 
negative photon momentum transfer squared, or photon virtuality, Q^, less than about 1000 GeV^. Above that value 
the QCD-Compton process becomes dominant, providing sensitivity to the strong coupling without being correlated 
to the gluon. Therefore, one can expect that the addition of jet measurements to the hitherto mainly used inclusive 
measurements alone will allow a reduction of the correlation between the gluon density and the strong coupUng, 
thereby making it possible to extract simultaneously both PDFs and the strong coupling from DIS data. 

HERAPDFS 



The idea behind the HERAPDFs is to use only HERA NC and CC data in the fits. Combined inclusive cross sections 
from the HI and ZEUS experiments are used, which have increased precision due to statistics but also due to 
intercalibration effects made possible by the two complementary detectors. Thus, precise combined data sets with 
total uncertainties between 1 — 2% (for NC) over most of the phase space are used. The systematic correlated and 
uncorrected errors are well controlled, such that the standard A^^ = 1 uncertainty criterion can be used. Since HERA 
data only include e^p data, there are no uncertainties due to for example deuterium or heavy target corrections. 

Two parametrizations have been studied: a 10 parameter fit which was used for the published HERAPDFl .0 [2] and 
for the prehtninary HERAPDFl. 5 [3], and a 14 parameter fit, for the preHminary HERAPDFl. 5f and HERAPDFl. 6 
[4]. The latter parametrization provides more flexibility, in particular for the gluon density at low x, which is preferred 
because of the additional jet data. At the Qq starting scale the PDFs are parametrized as follows: 

xg{x) = Agx'^!^ -{l-xfs -A'gX^'s -{l-xf's 

xuy{x) = A„„jt:^«" • (1 -x)*^"-' • (1 +D„^.jt: + £'„^.jt:^) and xd,.{x) = A,/,x^'''' • (1 - x)^*' 



xU{x) = Agx^iJ -{l-xfiJ and xD{x) = Ag.x^o • (1 -x)'^o 

In the 14 parameter fit Ag, A„,, and A^/,, are constrained by the momentum and quark sum rules. Furthermore it is 
assumed that Bq = Bjj, C'g = 25 and Aq — Ag(l — fs), where /j is the strange sea fraction at Qq. These choices define 
the central fit, as adding more parameters does not significantly reduce the x^- The parametrization uncertainty is 
estimated by adding one additional parameter at a time, changing the parametrizations and varying the starting scale 
2q around the default value of 1.9 GeV^. A model uncertainty is determined by varying the heavy quark masses, fs 
and Q^^„ around their default values. Q^^y^ specifies the minimal for a data point to be included in the fit. 

The fit of the PDFs is performed using the NLO DGLAP equations, employing the QCDNUM [5] evolution code. 
Charm and beauty quarks are treated as massive in the Thorne-Roberts variable-flavor-number scheme [6, 7]. The jet 
production cross sections are calculated with the fastNLO code [8], which efficiently convolutes the matrix elements, 
calculated using the NLOJet-H- program [9], with the fitted PDFs. 

Here is not the place to discuss this, but it should be mentioned that HERAPDFl .0 and HERAPDFl .5 fits have been 
also performed using evolution equations in next-to-next-to-leading order (NNLO). For the reader who wants to know 
more details and updates on the HERAPDFS, they can be found here [10]. 



DATA 

The data used in the extraction of the PDFs and the strong coupling as{Mz) are the combined NC and CC e^p 
and e^p inclusive cross sections from the HERA-1 [2, 11] and HERA-11 running periods [12, 11]. In addition, the 
following jet measurements from the HI and ZEUS collaborations were used: normalized inclusive jet cross sections 
for 150 < < 15000 GeV^ [13], corresponding to an integrated luminosity of 395 pb^', and inclusive jet cross 
sections for 5 < < 100 GeV^ [14], corresponding to an integrated luminosity of 43.5 pb^', and inclusive jet cross 
sections from two running periods in HERA-1 for > 125 GeV^ [15], corresponding to an integrated luminosity of 
38.6 pb^', and for 125 < < 10000 GeV^ [16], corresponding to an integrated luminosity of 82 pb^'. 



FIT RESULTS 



The PDFs as a function of x at ^ iq GeV^, when going from HERAPDFl. 5 (10 parameters) to HERAPDF1.5f (14 
parameters) with fixed as{Mz) = 0. 1 176, show an overall similar total uncertainty. The latter yields a softer sea density 
at high X and a suppressed d-valence density at low x. Keeping as{Mz) fixed and including the jet measurements, which 
leads to the HERAPDFl. 6 fit, one observes little difference in the PDFs and their uncertainties, only the high.ic gluon 
uncertainty is somewhat reduced. The respective figures can be found in [4]. More interesting results are obtained 
when the fits for HERAPDFl. 5f (no jets) and HERAPDFl. 6 (with jets) are performed with as{Mz) treated as a free 
fit parameter. The corresponding PDFs and their uncertainty bands are shown in Fig.l. Adding the jet data leads to a 
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FIGURE 1. (a) HERAPDFl .5f without and (b) HERAPDFl .6 with jet data included in the fit, as a function of x for = 10 GeV^ 
with as{Mz) as an additional free parameter in the fit, and (c) Ax^ for both fits. 

dramatic decrease of the low x gluon uncertainty, not only of the experimental uncertainty, but also of the model and 
parametrization ones. In Fig.lc the Ax^ distributions are shown as a function of as{Mz)- The HERAPDFl. 5f fit shows 



a shallow minimum, while the HERAPDF1.6 fit provides a strong constraint on the value of as{Mz)- The addition of 
the jet data successfully reduces the correlation of the strong coupling and the gluon density, for reasons discussed in 
the introduction. The increase in the uncertainty of the gluon density due to as{Mz) being a free parameter is not large 
when jet data are included. The fit yields the following value for the strong coupling: 

as{Mz) = 0.1202 ±0.0013 (exp.) ± 0.0007 (model/param.) 
±0.0012 (hadronization)+JJU°^^ (scale) 

The hadronization uncertainty is due to the corresponding correc- 
tions of the NLO jet cross sections. The uncertainty due to the scale 
is estimated by a variation of the renormalization and factorization 
scale by the conventional factor of 1/2 and 2 when calculating the 
jet cross sections. For the inclusive cross sections this uncertainty 
was assumed to be negligible. 

In Fig.2 the results on as (Mz) are displayed from HERAPDFl .6 
and as obtained from recent extractions by HI and ZEUS using 
only jet production data in DIS. The experimental and theoretical 
uncertainties are indicated separately. For HERAPDFl. 6 the un- 
certainty due to the PDFs is part of the experimental uncertainty, 
while for the jets only results it is part of the theory uncertainty. 
The values obtained using different methods and experiments are 
experimentally very precise. They are consistent with each other 
and with the world average. 

SUMMARY 

A NLO QCD fit with simultaneous determinations of the proton PDFs and as{Mz) was performed using HERA DIS 
data only. Combined HI and ZEUS inclusive NC and CC cross sections together with inclusive jet cross sections 
from the HI and ZEUS collaborations have been used. Including jet data in the fit and treating as{Mz) as a free 
parameter, dramatically reduces the correlations between the gluon density and the strong coupling compared to the fit 
without jets. The precision of the gluon density is improved, and an accurate and unbiased determination of as{Mz) 
is achieved, yielding a value consistent with the world average. 
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FIGURE 2. as[Mz) values from HERAPDFl. 6, 
from jet production in DIS, by the HI and 
ZEUS collaborations, and the world average. 



